Supported bilayers: combined specular and diffuse x-ray scattering 



o 
o 

(N 

-I— ' 
o 

O 



On 
O 



L. Malaquin^'^^ -p. Charitat^ and J. Daillant^ 
^Institut Charles Sadron, 
Universite de Strasbourg, CNRS UPR 22 
23 Rue du Loess, BP 84047 
67034 Strasbourg cedex 2, France 
CE A, IR AMIS, LIONS, CEA-Saclay, bat. 125, 
F-91191 Gif-sur-Yvette Cedex, France. 

A new method is proposed for the analysis of specular and off-specular reflectivity from supported 
lipid bilayers. Both thermal fluctuations and the "static" roughness induced by the substrate are 
carefully taken into account. Examples from supported bilayers and more complex systems com- 
prising a bilayer adsorbed or grafted on the substrate and another "floating" bilayer are given. The 
combined analysis of specular and off-specular reflectivity allows the precise determination of the 
structure of adsorbed and floating bilayers, their tension, bending rigidity and interaction potentials. 
We show that this new method gives a unique opportunity to investigate phenomena like profu- 
sion modes of adsorbed bilayers and opens the way to the investigation of more complex systems 
including different kinds of lipids, cholesterol or peptides. 



I. INTRODUCTION 

o , 

c/j , There is a great interest in determining the elastic properties and interactions between membranes at the 
microscopic scale. One attractive way to do it is to investigate their thermal fluctuations. This can be done for 
^ I example by recording the shape of vesicles (LUVs, GUVs) as a function of time and careful image analysis 
^ . The wavelength of fluctuations which can be investigated by this method is however limited to « Ifim. Shorter 
' ' wavelengths can be investigated by using x-rays. X-ray scattering has indeed been used to study unilamellar 
^ vesicles, and can be as precise as to demonstrate an asymmetry between the inner and outer leaflets Membrane 
Q ' fluctuations have however been mainly investigated in multilamellar vesicles [1, 0, The so-called Caille theory 
O . [1, 0] for x-ray diffraction from lyotropic liquid crystals enables in particular the determination of V kB where k is 
the bending rigidity of the membrane and B the stack compression modulus. Another advantage of the method 
Q is that the sample can be osmotically compressed in a systematic way. The structure of the sample can then 
^ ' be determined as a function of osmotic pressure 11, resulting in particular in n(a;) curves which can be fltted to 
' interaction potential models, usually taking into account hydration forces, van der Waals forces and the so-called 
Helfrich entropic interaction whose characteristic parameters can be evaluated. 

A signiflcant progress was the use of well orientated multilamellar stacks, either spin-coated onto substrates Q 
' I , or freely suspended where simultaneous fltting of specular and off-specular scattering allows an independent 
■ determination of k and B. However, all these methods intrinsically suffer from the large number of defects in the 
sample, leading for example to inconsistencies between the temperature dependence of the elastic parameters and 
the observation of an unbinding transition 



This issue can in principle be solved by looking at supported bilayers. When obtained by vesicle fusion it 
^ is likely that some defects will remain, but almost defect-free samples can in principle be obtained by using the 
Langmuir-Blodgett (LB) technique [l2l] . Their structure has been investigated in particular using atomic force 
,; microscopy [l3| and x-ray scattering. X-ray reflectivity has for example been used to determine the structure of 
d ' model lipid rafts in microfluidics cells p^ . and x-ray diffraction on lipid bilayers adsorbed on a substrate in water 
has recently been demonstrated p^ . 

One disadvantage of adsorbed bilayers is that strongly adsorbed s amp les would be very far from natural conditions. 
For this reason, polymeric cushions have been used as spacers [igI I17| |. Another possibility is to use a flrst bilayer as 
spacer, the second being less adsorbed on the substrate [l^ . A promising development is flnally the investigation 
of bilayer membranes spanning microfabricated holes [lij opening the possibility of investigating membranes in an 
asymmetric environment. 

In this paper, we consider adsorbed bilayers where the spacer is either a flrst lipid bilayer adsorbed on the substrate 
(double bilayers) [l^] or a mixed octadecytrichlorosilane (OTS) layer - lipid monolayer (OTS bilayers), where the 
OTS layer is chemically grafted on the substrate and the lipid monolayer is deposited using the LB technique p^ . 
A first study of such fioatin g b ilayers, mainly limited to the investigation of their structural and elastic properties 
has been recently published [20] and further extended to the determination of the interaction potential. The aim of 
this paper is to present the underlying theory leading to the full determination of structural and elastic properties of 
bilayer membranes and their interaction potentials. 
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We extend to more complex samples, in particular double bilayers, a first analysis of the height-height correlation 
functions of correlated interfaces by Swain et al. pll. [2^. Fluctuation spectra of free and supported membrane pairs 
have also been calculated in [13]. We show in particular that only a joint analysis of the specular and off-specular 
reflectivities allows a full characterization of the system. Specular reflectivity is more sensitive to the adsorbed bilayer 
whereas off-specular scattering is more sensitive to the floating bilayer. It consists in the systematic convolution of the 
substrate height-height correlations with susceptibility functions in order to appropriately propagate the correlations, 
and also takes into account thermal fluctuations. We then perform an efflcient numerical integration to obtain 
scattering cross-sections and intensities, taking precisely into account resolution functions. This method is quit e 
general and could be applied to similar problems like wetting films [23 | , polymer thin films or synthetic multilayers [25| | . 



II. SUPPORTED BILAYERS 



A. Free energy 




(b) (c) 

FIG. 1: Schematic view of supported bilayer systems on a rough substrate: (a) single supported bilayer; (b) mixed OTS- 
lipid supported bilayer; (c) supported double bilayer. We use the Monge representation to describe the membranes positions: 
r|| = {x,y) is the lateral coordinate, and z — Ui (ry) the position of the i*** interface (i=s (substrate), 1 (first layer), 2 (floating 
bilayer)). 



We consider a stack of almost flat membranes, with bending modulus Ki and surface tension 7^, supported on a 
rough surface. We use the Monge representation to describe the membranes as shown on Fig. [TJ r|| — {x,y), and z is 
the coordinate perpendicular to the substrate. We denote Ug (r||) the substrate position, and m (r||) the position of 
the i-th membrane (1 being the closest to the substrate). Each membrane interacts with the other components of the 
system through interaction potentials Us^i (r||) for the interaction with the substrate, and Uij (r||) for the interaction 
with another membrane. Following Canham 26] and Helfrich 27] we write the free energy of the system as: 

J^[u, (r||)]= Jd\ 

Considering small fluctuations of supported bilayers at equilibrium close to a substrate, we use a quadratic approxi- 
mation for the interaction potentials: 

Us,i (r||) = Us,i (ui (r||) - Us (r||)) ~ Ai (ui (r||) - Us (r||))^ 

U^.J (r||) = U,j (w, (r||) - (r||)) ~ S {u, (r||) - (r||))^ . (2) 

The various contributions to the interaction potential will be discussed in the next section. This approximation is no 
longer valid in the regime of large membrane fluctuations (see (28j). As usual for a system with harmonic coupling 
between degrees of freedom, it is worth rewriting the free energy of the system in Fourier space using: 
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I Q (t^V^ + K.A^) u^ (r||) + (r||)^ + U1.2 (r||) 



(1) 



zi,(r||)=5^zI,(g||)e''iii-"-|i, 

911 



(3) 
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where q|| = {q^, Qy). This transformation allows decoupling of the modes in the Fourier space leading to T = J^q^^ -^qw 
with: 

1 ^ 



where we use 5i(g||) = + 7ig| + Ki^jj^ with i — 1, 2. In the following, this formalism will be used to describe 
static and thermal deformations of membranes in different cases. 

B. Interaction potential 

The physical properties of lipid bilayers are the result of a competition between attractive and repulsive molecular 
interactions (van der Waals, electrostatic, hydration). The membranes are also subject to thermal fluctuations which 
lead to entropic repulsion [29l. fsol. [3l| . The van der Waals interaction U^dw^'^) between two membranes of thickness 
b separated by a distance z can be written: 

where the Hamaker constant H is on the order of fc^T and depends both on the lipids and on the solvent. At short 
length scales (less than 1 nm), bilayers separated by a distance d experience an exponentially decaying repulsive 
hydration force, the microscopic origin of which has been the matter of intense debate [H, [H, [H, [s^ . We write the 
hydration potential 

Uhyd = PhdhEXp (^~^^ ' (6) 

where dh — 0.3 nm ^ is the hydration length and 4 • 10^ Pa < < 4 • 10^ Pa is the hydration pressure. 
Finally, as first highligthed by Helfrich [365 > additionnal entropic contribution due to the confinment of the fluctu- 
ating membranes needs to be taken into account. Helfrich described the case of purely steric interactions ^36.] . leading 
to the well-known expression: 

This expression was extended in a phenomelogical way to tense bilayers (37| . Podgornik and Parsegian proposed a 
extension of this theory to soft-confinement, including direct interbilayer interactions [38| . 



III. CORRELATION FUNCTIONS 



In this section we calculate the height-height auto and cross-correlation functions of the membranes. These cor- 
relation functions will be used in Sect. IIVI to calculate the x-ray scattering cross-sections. As first described by 
Swain and Andelman '22*1, a supported bilayer deposited on a rough substrate adapts its equilibrium shape to 
the substrate roughness, and for finite temperatures, thermal fluctuations come into play. We will start by describing 
the substrate static roughness before considering different experimental systems (Fig. [T}: single supported bilayer, 
supported bilayer on a mixed OTS-lipid layer and supported double bilayer. In all cases, we first give the static 
correlation functions, generalizing the Swain and Andelman approach, and then describe the thermal contribution to 
the correlation functions. 

A. Static correlation functions 

1. Silicon substrate 

Like a wide range of rough surfaces, our silicon substrates can be described using self-afline correlation functions 

{uMus{rn))=<jle-(^^y"\ (8) 
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as, is the roughness ampHtude, < ff^ < 1 is the roughness exponent which describes the overall shape of the 
correlation function and is the characteristic length of the roughness (for example, distance between scratches...). 
For distances smaller than ^s, the correlation function goes to al as cr^(l — (.sY^') and decays to as a stretched 
or compressed exponential for distances larger than ^s. 



2. Single supported bilayer 

We first consider an almost flat single bilayer, with a bending modulus ki and a surface tension 71, interacting with 
a rough surface (Fig. [Ua)). The free energy is given by Eq. ([4]) with N=l. To determine the equilibrium state, this 
equation has to be minimized against ui{q)., yielding ui{q) = {Ai/ai{q)) Us{q)- The roughness spectrum is given by 

Hi HI: 

{Mq\\)M-q\\)) - -#t(I^s('?ii)P> = ^ -^(l^.(<zii)P> (9) 



^1 +7i9f| + t^iQii 



ai[q\\) /li + 7i(7|| + Ki^ii 

For low bending rigidity and/or tension, the membrane is soft enough to follow the substrate roughness and we have 
(|ui(q||)p) ~ (|us((7||)p). On the opposite, if the rigidity or the membrane tension are large, (|{ti(g||)p) — > as the 
membrane is too stretched or rigid to follow the substrate. Two cases can be considered. If A = 7^ — AA^ki > 0, the 
membrane follows the substrate for wave- vectors qn < yMi/71, there is then a tension dominated regime in for 

\/ Ai/"fi < (7|| < ^71 and finally a bending rigidity dominated regime in I/qI for q\\ > ^^71/^1. If A < 0, we 
directly go from the potential dominated regime where the membrane follows the substrate to the bending rigidity 
dominated regime. 

In the real space, Eqs. ^ and (fTO| become convolution products: 

(ui(OK(r||)) - 5i(r||)(^(M,(0K(r||)) (11) 
(ui(O)wi(rii)) = /ii(r||)® (M,(0)u,(r||)), (12) 

where the susceptibilities gi(?'||) and /ii(?'||) describe how the membrane self- and cross-correlation functions correlate 
to the substrate heght-height correlation function. 5i(r||) is the inverse Fourier transform of A^/a^(g||) and ^i(?'||) the 
inverse Fourier transform of Ai/di(q^\). In order to calculate gi{r\\) and hi{r\\) we need the roots of ai(9||) given in 
appendix I VII B D in order to decompose gi and hi in partial fractions. In the following we denote them as gin and q2\\. 
The Fourier transforms of polynomial fractions such as l/(?^ + ^In) or l/(gj| -|- 9|||)^ are proportional to the modified 
Bessel functions of the second kind KQ{qj^^r^^) and Ki{qj^^r^^) (see appendix IVII A[) and we obtain: 



2 2 

XiKo (qiiirii) + X2K0 (gzirii) + i^i (gi||r||) + ^-^i^i (qaFn) 

^ 9111 ^ 9211 



(13) 



^i(^ll) = [mKa {qi\\r\\) +V2K0 {q2\n)] , (14) 

where the rji and Ai are given in App. IVII Al 



3. Two membranes near a wall 

In this section, we treat the case of two supported membranes (N=2, Fig. [Ijc)). The minimization of the free 
energy against ui{q\\) and U2{q\\) gives: 

~ f ^ Aid2{q\i) + B (Al + A2) 

ai(q\\)a2{q) + B [ai[q\\) + a2{q\\)) 
~ I ^ A2ai(g||) + S(Ai+A2) ^ 

"2(911) = ~ I ^~ I \^-a(~ ( ( (16) 

ai (9|| )a2 (9|| ) + S (ai {q\\ ) + 02 (g|| ) j 
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Five correlation functions are now needed to describe the membranes: (i) the self height-height correlations of each 
membrane (ui(0)ui(r||)^ and (^2(0)^^2(^11)); (ii) the height-height cross-correlations between the membranes and the 
substrate (us(0)Mi(r||)) and (^Us{0)u2{r\\)); (iii) and the cross-correlation between the two membranes (ui(0)u2(r||)). 
These functions can be expressed as a convolution product of Eq. ([H]) with a polynomial fraction in . Applying the 
method of the previous section, we first determine the roots of the polynomial equation of degree 8, ai((7||)a2((j'||) -I- 
B (51(911) + 02(911)) — and perform the partial fraction decomposition. Like in the previous section, the correlation 
functions can be expressed as combinations of the Bessel functions Kq and Ki: 



(ui(0)M2(r||)) = fu2 {r) {us{0)us{^)) 
(Mj(0)w,(r||)) = h,^, (r) (g) {us{0)u,{r\\)) , 



(17) 
(18) 
(19) 



with: 



/i^2 (r||) 

5*(ni) 

hi,s {r\\) 



1 A1A2 
27r K1K2 



E 



2- U 



TjKo{qj\m] 



27r K, 



1 Ki{qj\m) 



1 2 



Kiiqjin) 



(20) 



(21) 



(22) 



Eqs. (pij) and ([^ are the generahzation of the Eqs. and (fH)l for double bilayers (see App. IVII B 21 for the 

expression of the various coefficients). Whatever the number of membranes, these correlation functions will always 
have the same shape. It is worth noting that for K2 = 72 = ^2 = S = 0, we correctly recover the correlation function 
of the single bilayer obtained in the previous section. 

In the specific case of symmetrical bilayers (71 = 72, ki = K2 and Ai — A2), we obtain 9311 — qi^^ and g4|| = g^n, and 
the parameters Xi_j rjij, Tj and i^j can be easily calculated: 



1^1 



T2 = Ai. 



A, 



'2|| 



^l|i 



and A1.2 = A2,2 = -n 



and 



1^2 = vli = vli 



Vl,2 = V2,2 = ^i- 



(23) 



The remaining terms, A^^s, Ai,4, ry^^s, 77^4, 1^3, v^, and are all equal to zero. 



B. Thermal fluctuations 



In addition to the static roughness induced by the substrate described in the previous section, the supported 
membranes undergo thermal fluctuations which we discuss in this section. The total displacement of a membrane 
is Ui{q\\) = Ui^s{(l\\) + M'i,t/i(9||), where Ui^s{q\\) represents the static position of the i-th membrane, and Ui^th{(l\\) its 
thermal fluctuations. Assuming that there is no correlation between {ii_s(9||) and Mj,t/i('?||) we can separate the free 
energy into a static part J^s,q^\ and a thermal one Tth,q^\ ■ 



1. Single bilayer 

For a single bilayer, the thermal part of the free energy is J-th,q\\ = ^Oi (Q|| )\uth,i{(i\\ ) By applying the equipartition 
theorem one obtains: 

{uth,i{q\\)uth,i{'q\^)) = = -. — , ^^2 , 4- (24) 

ai{q\\) Ai+7ig^ + Kig* 
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In the real space, the height- height correlation function is simply given by (20|: 

1 ks". 
27r Ki 



1 kj^T 

{uth,ii^)'^th,i{r\\)) = — [ai/-i:o('?i||f||) + a2-?^o('Z2||''||)] ■ (25) 



Coefficients ai and a2 are given in App. I VII D II 

In the case of a mixed OTS-lipid double bilayer, we only have to replace ai(5||) with 6(q||) in Tth.q^^ 



2. Double bilayers 
The free energy Eq.(IT|) for a double bilayer can be written: 

1 ^ 

•^t'i,9ii = ■7^'^{^M\\) + B) \uthA(l\\)? ~ ButhAq\\)uthA~(l\\)- (26) 



1=1 



The fluctuation modes of the two membranes uth,i{c[\\) and uth,2{(l\\) are coupled in Eq. ([26|l and the Hamiltonian of 
the system, 

~ai{q\\)+B -B/2 
-B/2 a2{qi\) + B 

needs first to be diagonalized in order to apply the equipartition theorem. The calculation is detailed in App. IVII CI 
where we find: 

/l~ f m2\ I, rj. a2{q\\)+B 



ai(gi[)a2(g||) + B (ai(gj|) + a2{q\\)) 

ai(g||)a2(g||) +B (ai(g||) + a2{q{)) 
(ut/i,i(9||) • ut/i,2(-9|i)) = ksT- 



aiiq\\)a2{q\\) + B [ai{q\\) + a2{q\\)) 

To determine the inverse Fourier transform of Eqs. (|27p . one first needs to determine the roots of the equation 
of degree 4, ai{q\\) + B — 0, and the roots qq of the equation of degree 8, ai(g||)a2(g||) + B (ai(g||) + a2{q\\)) = to 
perform the partial fraction decomposition. It is worth noting that the equation of degree 8 is the same as in the 
static case. Like before, the correlation functions then can be expressed as a combination of the Kq and Ki modified 
Bessel functions: 

{uthA^)uthAr\\)) = -^^y2'^^,jKaiqJ\\r\\) 
k T B ^ 

{uth,ii0)uth.2{r\\)) = -^:5-yi^i^o(gjini)' (28) 

K1K2 

J = l 

where the parameters are given in App. IVIID 21 

In the particular case of uncoupled bilayers (B — 0), these equations can be simplified, leading to: 

9i|| = 93|| = = f32,i; 9211 = 9411 = 01,2 = (32,2 

1 

= a2 ,1 = 5-; q;i,2 — OL2,2 = — 

9211 - 9111 

All the remaining parameters a^^a, aj^4 and li are equal to zero and we recover the results of [20| : 

/b T 1 1 

(uth,i{^)uth,i{r\\)) = — — 2- (■^o('72|ri|) - ^''o('7i||''||)) 

Ki An Q2 II ^1 II 



{uth,i{0)uth,2ir\\)) = 0. 



7 





4.0 10"^° 




^ ^ «-20 

3.5 10 




3.0 10'^° 




2.5 10"^° 








2.0 10'^° 






<u(r 


1.5 10 




1.0 10'^° 




5.0 10"^^ 




0.0 10° 





4.0 10"^° 




3.5 10"^° 




3.0 10"^° 




2.5 10"^° 






o. 


2.0 10'^° 






<u(r 


1.5 10=^° 




1.0 10"^° 




5.0 10"^^ 




0.0 10° 



substrate 

membrane-membrane (static) 
substrate-membrane 




4.0 10 



10"' 



substrate 
7= 0.1 mN/m 
7 = 1 mN/m 
7=10 mN/m 



10" 




3.5 10 
3.0 10 
2.5 10 
2.0 10 
1.5 10 
1.0 10 
5.0 10 
0.0 10' 



20 



20 



20 



20 



20 



20 



20 



10"" 






4.0 10"^° 




3.5 10"^° 




3.0 10"^° 




2.5 10"^° 








2.0 10"^° 






V 


1.5 10"'° 




1 .0 1 0"'° 




5.0 10"'^ 




0.0 10° 



substrate 
U" = 10^J/m^ 
U" = 10 , J/m* 
U" = 10 "j/m^ 
^■ = 10^°^ 




FIG. 2; (a) Static correlation functions for a bilayer with 7=1 mN.m"^, k = 30 ksT, Ai = 10^^ J.m~* on a substrate with 
ffa = 0.19 nm, ~ 170 nm and Hs = 0.5. Influence of the different parameters on the bilayer auto-correlation function: (b) 
bending rigidity ki , (c) membrane tension 71 , (d) second derivative of the effective potential Ai . 



C. Discussion 



Fig. [5] sho-ws the effect of the elastic properties of a single supported bilayer on its coupling to a rough substrate 
("static" correlation function). The substrate-substrate, substrate- membrane and membrane-membrane correlation 
functions for the bilayer are plotted on Fig[2^. In our model, the correlation functions of the membrane (auto- 
correlation and cross-correlation) are al"ways smaller than that of the substrate and the cross-correlation function is 
always comprised in bet"ween the substrate and membrane correlation functions, "whatever the elastic parameters. In 
all cases, the bilayer foUo-ws the substrate at large lengthscales, "whereas the short lengthscale behavior depends on 
the potential and the elastic properties. 

Fig. [2}3 sho"ws the membrane-membrane correlation function for two relevant values of the bending rigidity, k = 30 
ksT (fluid phase) and k = 300 k^T (gel phase). The r.m.s. membrane roughness (limit of the correlation function 
"when r|| — > 0) decreases for increasing k as a more rigid membrane is less free to follow the substrate roughness. 
For the same reason, the cutoff in the susceptibilities 51 (^'n) and /ii(r||), Eq. p2p . increases with k (from 70 nm for 
K = 30 ksT, 7=1 mN/m and A = 10^^ J.m"" to 209 nm for n = 300 k^T). As noticed in Sect. IIII A 21 in between 
K = 30 ksT and k = 300 k^T we move from a regime where the cutoff Tc is determined by the tension and rigidity 
{vc — 27r-y/ k/7 when A = 7^ AkA is positive) to a regime where it is determined by the potential and rigidity 
(re = 27r(K/A)i/4 when A < 0). 

Similar trends are observed with increasing tension (Fig. [2fc) and potential (Fig. [2]i). For example, the cutoff is 
equal to 22 nm (roughness ) for 7 = 0.1 mN/m k = 30 k^T, and A = 10^^ J-m^**), 70 nm (roughness ) for 7 = 1 
mN/m and 118 nm (roughness ) for 7 = 10 mN/m, as a more stretched membrane cannot follow the substrate at 
short lengthscales. Here again, the cutoff is determined by potential and rigidity for 7 = 0.1 mN/m and by tension 
and rigidity when 7 = 1 mN/m or 7 = 10 mN/m. 

The thermal correlation function of a supported bilayer close to a substrate is plotted on FigO using the same set of 
parameters. The thermal fluctuations can be larger than the substrate roughness if the interaction potential is soft 
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FIG. 3: (a) Thermal correlation functions for a bilayer with 7=1 mN.m ^, k = 30 ksT, Ai = 10^^ J.m~4. (b) Effect of the 
bending rigidity k on the membrane correlation function, (c) Effect of the tension 7 on the membrane correlation function, (d) 
effect of the second derivative of the effective potential Ai on the membrane correlation function. 



enough, or if the membrane tension and rigidity are small enough. This is in particular the case in the fluid phase 
(Fig. [3^) , "where the thermal roughness is almost one order of magnitude larger than the static roughness for realistic 
parameters. The elastic properties of the membrane have the same effect on static and thermal correlation functions. 
The average roughness of the bilayer and the cutoff increase -when the bending rigidity (Fig. the tension (Fig. 
[5]:) or the second derivative of the effective potential (Fig. decrease. When A = 7^ — AkA > 0, the membrane 
tension comes into play, and the correlation function has the characteristic logarithmic decay due to its long-range 
effects. 



IV. SPECULAR AND OFF-SPECULAR INTENSITY SCATTERED BY A MEMBRANE 



Specular and off-specular X-ray reflectivity can be used in order to determine the membrane structure and properties 
at submicronic length scales. Specular reflectivity allo"ws one to determine the structure of the sample perpendicular 
to the average membrane plane whereas off-specular reflectivity also gives access to the elastic properties of the mem- 
brane (correlation functions, tension, bending energy) and to the interaction potential. 

In this section, we show how the scattered intensity can be derived by using the correlation functions obtained in the 
previous section. The geometry of the experiment is defined on Fig. [J] The grazing angle of incidence is 0in- Osc is 
the angle of the scattered x-rays in the plane of incidence (jp normal to it). 

In specular refiectivity experiments, the sample is rocked around the specular condition for every detector position 
(thus keeping the normal wave-vector transfer Qz approximately constant) in order to record and subtract the back- 
ground. In off-specular reflectivity, the grazing angle of incidence is kept fixed {Oin = 0.7 mrad) below the critical angle 
for total external reflection 9c « 0.85 mrad for the silicon-water interface at 27 keV, as this allows easy background 
subtraction (20| . and 9sc is scanned in the plane of incidence. 

In order to calculate the intensity scattered by the sample, we first calculate the differential scattering cross-section 
da/dfl which is the power scattered by unit solid angle Qd in the direction ksc per unit incident flux in the direction 



qx 

(a) 



(b) 



FIG. 4: (a) Schematic view of tlie experimental setup for specular and off-specular reflectivity and (b) Fourier space trajectories 
in specular and off-specular reflectivity experiment. The grazing angle of incidence is Bin. 9sc is the angle of scattered x-rays 
in the plane of incidence (t/i normal to it). 



ki„ (Fig. 2^). The scattered intensity is then calculated by integrating over the detector sohd angle, 

d(7 



I in) 



(29) 



Differential scattering cross-sections are large close to the critical angle for total external reflection and multiple 
scattering cannot be neglected. In other words, the simple kinematic Born approximation is no longer valid and 
a better approximation must be used. Here we use a simplified "Distorted Wave Born Approximation" (DWBA) 
[20I [25l . Ho], which is a perturbation theory using as reference state a perfectly flat silicon/water interface. As the 
electron density of lipids is close to that of water, this is in fact an excellent approximation. 



A. Specular reflectivity 



Within this approximation, the specular reflectivity can be written [20|, |25| : 



R{q,) = RF [qz) 



5p{z) 



PSi — PH2O 



Rf {qz) 



1 



PSi — PH2O 



^ 1 e^^-dz 
dz 



(30) 



where p is the electron density and dp is the difference in electron density between the real system and the reference 
state. 



B. Off-specular reflectivity 



1. Differential scattering cross-section 



Within the DWBA the perturbation part of the differential scattering cross-section is given by [20| : 

2\ 



J dr||e'i '-|i<5/9(z) 



(31) 



where — 2.8 x 10 -^^ m is the classical radius of the electron, <™ — tH20,Si{^in) and t'^'^ — tH2O,Si[0sc) are the 
Fresnel transmission coefficients for the silicon/ water interface, i*" is a good approximation to the incident field 
scattered by the interface, and f"^ describes how this field propagates to the detector. (ei„ • egc)^ is the polarization 
factor. In our case (ei„ • Gsc)^ — 1 [40j . As our sample is composed of a substrate and of two bilayers, we can split 
the electron density in three terms: SpsxCo for the substrate and (5pMi and (5pM2 for the two membranes. Each of these 
terms depends on rii because of the static roughness and of the fluctuations. Wc have 



5p {z, r||) = Spsnh [z - Us) + SpMl {Z - Ui,th - Ml,st) + 5pM2 {z - U2,t/i - U2.st) , 



(32) 
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where Us, Ui^st and Ui^th depend on r||. Jpsub = Psub - PH2O for < z < if > and Spsuh = Ph^O - Psub for 
Ms < 2; < if Ms < 0. Spuii — PMi ~ PH2O with i — 1,2 are the electron density profiles of membranes relative to water. 
After change of variables z — Ui^th — Ui^st z, we obtain: 



dr\\ 



(33) 



where Spsuh (qz), Spui (qz) and SpM2 (qz) are the substrate and bilayer form factors. Assuming that there is no cross 
correlation between the thermal and static fluctuations of the membranes and between the thermal fluctuations of 
the membranes and the substrate roughness, we get: 



{[ui,th (0) + u^^st (0)] [ui^th (r||) + Ui^st (r||)]) = (ui.th (0) Ui^th (r||)) + (ui^st (0) Ui^st 



(34) 



We then proceed by developping Eq. ([55]) . Using the classical result (e*^^) — e * for a Gaussian variable z, we 

can rewrite the differential scattering cross-section as: 



+ \SpMi ((j^)|^e~''-('^i-"'+''i'=')e'''("i'"'^°^"i-"' (""11 ))e'''("i'=''°^"i'=' ('■")) 

+ \SpM2 (g^)|^e"''-(''''"•+'^2,=t)e9^("2."'(0)"2,t^^(r||))g9"(«2,et(0)«2, = t(r||)) 

+ iSplii (qz) SpM2 (qz) + SpMl (qz) Sp*M, (g.)) e"^?^ (-1-+-L,*+-L.+-L0 ("^'^ ("^"M. (r|| )) ^9= («2,.* (0)«i,.. (r,, )) 

+ iSp*Ml (Qz) <5ftub (qz) + SpMl (qz) 5pU (9.)) g-^^^' K.-+<-+-ne'^' W"^^- (^1 )> 

+ iSp*M2 (qz) Spsub [q.) + 6pM2 {qz) SpU (qz)) e-^«^'(-^.-+-^^-+-?)e«^'("=(°)"^-'("-ii)>l 



(35) 



where crf ,,^ = ^Ui^st (0)^y and 

'^i,th — i^i.th (0) I are respectively the static and the thermal siuface roughness of 

the i-th membrane. A is the illuminated area on the sample [41,]. Finally, rewriting — 1 + — 1, specular and 
off-specular parts of the perturbation in the scattering cross-section can be separated: 



g(g,)(27r)2 5(q||) + / dr||e^'i|l"^ii/ (r,,, 9.) 



(36) 



with: 



-t- 
-I- 



/(r||,g,) = |p3„b(g.)re-^^^'^^(e'^^(-(°)-(-^"))- 

\PM2 (g^)|^e"'^'('^'-*''+''^-') |^e'''("2,t;.(0)u2,th(r||))gg^(«2, = t(0)«2,et(r|i)) „ 
+ {PMl{qz)pM2{qz)+PMl{qz)plA2{qz))e-"^'i"^<^-+"^^^^^^ (^e«^ («2..(0W,, (r,, )) («2,.* (0)«,,.. (r„ )) 

+ (PMl {qz) Psub {qz) + PMI {qz) pU {qz)) (e«.?('-(0)"i,-(r,l)) _ 1) 

+ (PM2fe)Psub(g.) + PM2(g.)p:ubfe))e-^'^('^--+<-+^^Ue'^(^^^ 



The g{qz) function will not be used in the following as we prefer to use Eq. pO|) which directly gives access to the 
total reflectivity and not only to the perturbation part. 
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2. Off-specular intensity 



Resolution effects can be taken into account by using a resolution function TZes {9sc,'4'), which is equal to one if 
'sc: ^) points to the detector area and outside, 



/ (Si) M + ( „ S 

* * J \ / spec \ / off— spec 



(38) 



hi and Wi are the incident beam height and width so that lo/hiWi is the incident flux. Contrary to the specular 
case where resolution effects amount to a convolution, the off-specular intensity is proportional to the detector solid 
angle. The experimental resolution is a unit rectangular function centered on the detector, of width in the plane 
of incidence and Aip perpendicular to it. It is in fact easier to proceed in the Fourier space where we have: 



qx = ko (cos6li„ - cos 6'sc) dqx = ko sm6scd9 

Qy = ^0 sinV' — fcoV' =^ dqy = k^dip. 

For off-specular scattering, Eq. (j38p becomes: 

^0 



(39) 
(40) 



^off— spec 



hiWikQ sin 



da\ 



dfl J CI 

/ off— spec 

iWi k^smOsc 



TZes ( q|| ) dq| 



dxjdyj dq, J dqye'''^^e'''yyf (r,, , q,) Ties (q,,) . 



Since, under our experimental conditions, the slits are widely open in the y direction, we have TZes{qp 
The integration over qy and y is performed using J dqyC^'^^^ — 2tt5 [y). We first get: 



-^off — spec 



MWi fcp sin 



dx / dg,e'?^V((a;,0),(7,)7^es(<7,). 



(41) 
TZes{q^). 

(42) 



We further proceed by approximating the resolution function along qx as a Gaussian function of width Aqx — 
kasinOgcAOsc centered in q^o'- 



Ties {qx) 



The term 1/\/2tt is used to normalize the function. One obtains: 



-^off — spec — 

Integrating over qx gives: 



2iWi k^sinOsc 



dqx / da;e'«-o^e*9-^e ^^/(x,g,). 



(43) 



(44) 



(45) 



(46) 



which is the main result of this paper. An efficient method for the numerical integration of Eq. (j46[) is given in App. 

lylTEl 



dqxC 



We finally obtain: 



2 „ + oo 



V. EXAMPLES 
A. Experimental 



In this section we discuss several examples where the results from the previous sections are used to fit experimental 
data. All the experiments reported here used a 27 keV x-ray beam (wavelength A = 0.0459 nm) at the CRG-IF 
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beamline of the European Synchrotron Radiation Facihty (ESRF). The scattering geometry is described in figure 
5) The monochromatic incident beam was first extracted from the polychromatic beam using a two-crystal Si(lll) 
monochromator. Higher harmonics were eliminated using a W coated glass mirror, also used for focusing. In all 
experiments, the incident beam was 500 /im xl8 /im (W x H). The reflected intensity was defined using a 20 mm 
x200 /im (W X H) at 210 mm from the sample and a 20 mm x200 /im (W x H) at 815 mm from the sample and 
recorded using a Nal(Tl) scintillator. 

Bilayers and double bilayers of L—a l,2-distearoyl-sn-glycero-3-phosphocholine (di-Cis-PC or DPSC) from Avanti 
Polar Lipids (Lancaster, Alabama, USA) were prepared by first depositing a bilayer by two classical Langmuir- 
Blodgett (LB) depositions (vertical sample). We used super-polished (< lA roughness) silicon substrates (SESO, 
Aix-en Provence, France) of surface 5x5 cm^ and 1 cm thick (to ensure planarity). The floating bilayer was then 
prepared by a LB deposition, followed by a Langmuir-Schaeffcr (LS) deposition (horizontal sample) [l2|- The transfer 
rates are measured with 0.01 precision and 0.02 statistical dispersion. 

Mixed OTS-lipid bilayers were prepared by first coating the substrate with an octadecyltrichlorosilane (OTS) as 
described in [i9|. Then, 2 monolayers of DSPC are deposited by Langmuir-Blodgett technique and a last one by 
Langmuir-Schaeffer technique. 

The samples are then inserted into a PTFE sample cell with 50^m thick windows which is put in an aluminum 
box and thermalized using a water circulation bath. Sample are heated by steps with a feedback on the temperature 
measured inside the sample cell using a PtlOO resistance. Specular reflectivity is obtained by rocking the sample for 
each angle of incidence (q^ scans for approximately constant ) of incidence in order to subtract the background. 
Off-specular reflectivity is measured at a constant grazing angle of incidence of 0.4 mrad bellow the critical angle to 
total external reflection at the Si- water interface (0.7 mrad). Background subtraction is described in detail in [2Q] . 

B. Results and discussion 

We first determined the silicon substrate correlation function parameters, Eq. ([5]), by fitting the calculated specular 
and off-specular intensities to the experimental data for the substrate- water interface. Fig. [5] shows the effect of these 
parameters. As expected, the off-specular scattering is more sensitive to the shape of the height-height correlation 
function than the specular reflectivity, which is in turn more sensitive to the electron density profile (thickness and 
density of the Si02 layer). The best set of parameters is (Ts = 2 ± 0.2A, = 1 ± O.l/im and Hs — 0.5 ± 0.05, in 
good agreement with AFM characterization. Whereas a change in the roughness as results in an overall shift in the 
scattered intensity in logarithmic scale, a change in the cutoff results in a shift in the crossover from a constant 
intensity regime at low q^ to a l/q^.'^'^^'' regime at larger q^ values. 

The experimental results and best fits for off-specular and specular reflectivity are given in Figs. [n]and[7]for a single 
bilayer at T = 20°C, a double bilayer at T = 49°C and T = 62°C and a mixed OTS bilayer at T = 42.9°C. Good 
flts could be obtained in all cases for both specular reflectivity and diffuse scattering. In all cases, the off-specular 
scattering from the bare substrate provides a baseline for the contributions from the bilayers (Fig. [7]). Electron 
densities which are represented on Fig. [5] were modelled using the so-called 1-G gaussian model consisting in one 
gaussian for the head groups, a fiat part describing the chains and a Gaussian methyl trough [13, Parameters are 
given in App. IVIIFI 

Whereas the specular curves (Fig. [S]) are mainly sensitive to the average structure of the membrane, the off-specular 
signal (Fig. [7]) is both sensitive to the elastic properties of the bilayer and interaction potentials, and to the structural 
parameters (because both q^ and qz are varied in a detector scan. Fig. [ijb)). In particular, whereas the shape of the 
specular signal results from interferences between the qx ^ components of the electron densities where all interfaces 
are fully correlated, the diffuse signal is directly determined by the strength of the correlations between interfaces at 
a given q^, Eqs. and (|57|) . leading to a reduced contrast. This is demonstrated on Fig. [51 where in addition to 
the reduced contrast, a contrast inversion between specular and diffuse refiectivities (maxima in the specular intensity 
correspond to minima in the diffuse scattering) can be seen at low q^ as bilayer-bilayer correlations dominate over 
substrate-bilayer correlations. As a consequence, the substrate is not seen in low q^ diffuse scattering, leading to the 
contrast inversion. 

It should also be noted that the experimentally measured refiectivity signal consists of a truly specular signal and 
a non-negligible diffuse contribution, both shown on Fig. [SI as the detector slits have a finite opening. The latter is 
relatively more important for large qz values as the specular signal decays as 1/g^ exp(— g^cr^), Eqs. (j30p . where a is 

a typical roughness value, whereas the off-specular signal decays as l/g^"*"^-^" ~ l/ql/'^^'^% Eqs. (|^ and (|57)l . 
Due to the exp(— g^cr^) decay in the specular signal, the contribution from the more strongly fluctuating floating bilayer 
to the specular intensity vanishes at large qz- This is in contrast to the diffuse scattering where larger fluctuations 
make a larger contribution. Accordingly the more strongly fluctuating floating bilayer makes the largest contribution 
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10^ 10^ 



FIG. 5: Top: effect of the average roughness (Ts on the off-specular intensity. Middle: effect of the correlation length ■ 
Bottom: effect of the roughness exponent Hs- Inset: best fit of the specular reflectivity. 



at large > 10 ""m^^ (see in particular the double bilayer in Fig. [7]), whereas the adsorbed bilayer and the substrate 
- adsobed bilayer cross-terms make a larger contribution at lower . Combining specular (Fig. ^ and off-specular 
(Fig. [7]) reflectivity is therefore essential in order to obtain a good sensitivity to the structural parameters of both 
bilayers (Fig. [8]). 

The sensitivity of the diffuse scattering to the tension, bending rigidity and second derivative of the potential is 
demonstrated on FigdOlfor the mixed OTS-lipid sample. Unless otherwize specified, on FigfTOl k = 300 ks T which 
is a typical value for a lipid bilayer in the gel phase, 7 — 0.4 mN/m, U" — lO^^ '^^ J.m^^ and IH2O ~ 2.3 nm which 
provide the best fit to the experimental curve. A key region in the scattering curves is the fringe close to (7|| w 5 x 10^ 
m~^ or Qz « 1.2 X 10^ m~^, see also Fig. [9l which results from the interference between x-rays scattered by the two 
bilayers. Its location is therefore extremely sensistive to the water layer thickness (Fig. [TUH) and its contrast is most 
sensitive to the second derivative of the potential (Fig. [TUb). Too large a tension kills this correlation between the 
two bilayers. However, the main effect of a smaller (resp. larger) tension is to increase (decrease) the scattering 
at low g|| and decrease (increase) it at large due to the larger (smaller) roughness (Fig. [TUk). Finally, the main 
effect of the bending rigidity is at large (7|| (Fig. [TUb). As all these parameters affect different regions of the curve in 
different ways, they are accurately determined in the fits without too much coupling. 
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FIG. 6: Specular reflectivity and best flt (full line) for a supported bilayer at 20° C (a), a double bilayer at 49° C in the gel phase 
(b), 62° C in the fluid phase (c) and a mixed OTS-lipid bilayer at 42.9° C in the gel phase (d). The truely specular contribution 
is indicated using dashed lines and the contribution of the diffuse scattering in the specular direction is indicated using dotted 
line. 



Regarding the single bilayer, the best fit yielded [/" = ioii-2±0-5 j.ni"^ for the second derivative of the 
potential and a membrane tension 7 « 70 ± 20 mN.m^^. With this large 7 value, the fit is insensitive to the 
bending rigidity k. This large tension is much larger than lysis tension (w 5 — 20 mN.m~^) [44] and can be 
interpreted as an effective values describing lipid profusions rather than a usual tension. Profusions are local 
independent motions of lipids, which should be the only possible motion for this strongly adsorbed membrane, 
they tend to increase the local surface area and can therefore be described using an effective microscopic tension 
as suggested by Lipowsky and Grothehans [isj. Lindahl et al. [1^ simulated the behavior of a fioating bilayer 
found 7p — 50 mN.m^^ which is in the same order of magnitude as the surface tension we determine in our experiment. 
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FIG. 7: Diffuse scattering and best fit (full line) for a supported bilayer at 20°C (a), a double bilayer at 49°C in the gel phase 
(b), 62°C in the fluid phase (c) and a mixed OTS-lipid bilayer at 42.9°C in the gel phase (d). The contribution of the first 
bilayer close to the substrate is indicated using dashed lines and the contribution of the floating bilayer is indicated using dotted 
line. 



Regarding double bilayers and mixed OTS bilayers, we obtain k w 200 — 300 kg T in the gel phase which is in 
good agreement with the results found by Pecreaux et al on giant unilamellar vesicles [ij and less than 50fcBT in 
the fluid phase. Interestingly, we consistently obtain the efl^ective tension of the first bilayer close to the substrate 
to be one order of magnitude larger than that of the floating bilayer. Whereas this can be attributed to profusion 
modes for double bilayers, a large effective tension is also not surprising for a mixed OTS-lipid bilayer where the 
OTS bilayer is polymerized and anchored on the substrate. The tension of the floating membrane, 7 < 2 mN/m, 
is larger than for vesicles but consistent with the fact that the membrane is still supported even if we consider it 
as a "floating" bilayer. Regarding the second derivatives of the interaction potential, it is interesting to note that 
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Distance [A] 

FIG. 8: Electron-density profile of a di-Ci8-PC obtained by combined specular and ofT-specular reflectivity using a Gaussian 
model comprising one Gaussian for the head groups and one Gaussian trough for the methyl groups [4^. l43|. (a) Bilayer, (b) 
double bilayer 49°C, (c) double bilayer 62°C, (c) mixed OTS-lipid bilayer 42.9°C. 



the smallest coupling between the first bilayer and the substrate is obtained for the mixed OTS-lipid bilayer, which 
is consistent with a polymerized OTS layer grafted at a relatively small number of points. On the other hand, the 
largest coupling is obtained for the adsorbed bilayer in the fluid phase. The second derivative of the interaction 
potential between the floating bilayer and the substrate is on the order of (1 — 3) x 10^° J.m"'* in the gel phase for 
either the double bilayer or the mixed OTS-lipid sample and decreases by more than one order of magnitude in the 
fluid phase. Finally, the largest interactions are between bilayers. 



17 




FIG. 9: Specular reflectivity and diffuse scattering from a double bilayer at 49°C. x I/Io as a function of q^. The maxima 
in the specular intensity correspond to minima in the diffuse scattering at low qz. 



VI. CONCLUDING REMARKS 

The formalism developed in this paper provides a rigourous frame for analyzing scattering data from supported 
bilayers. Coupling of specular and off-specular scattering is essential for determining the structure of adsorbed and 
floating bilayers as well as their tension, rigidity and interaction potentials. Several systems have been investigated, 
absorbed bilayers, double bilayers, mixed OTS-lipid bilayers. Our data show that this new method gives a unique 
opportunity to investigate phenomena like profusion modes of adsorbed bilayers and opens the way to the investigation 
of more complex systems like charged membranes or more realistic models including different kinds of lipids, cholesterol 
[131 or peptides. 

More generally, our method based on a careful propagation of correlation functions provides an efficient scheme for 
tackling different systems like wetting films, polymer layers... where interaction potentials and/or elastic parameters 
could be accurately determined. 

Ackonowldgements: we wish to thank J.-S. Micha for assistance during the experiments, and G. Fragneto, F. Graner 
and S. Lecuyer for discussions. 
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FIG. 10: Effect of the elastic parameters of the floating bilayer, second derivative of the potential and water layer thickness on 
the diffuse scattering from the mixed OTS-lipid bilayer at 42.9°C. (a) Effect of the surface tension of the floating bilayer. (b) 
Effect of the bending rigidity, (c) Effect of the second derivative of the potential, (d) Effect of the thickness of the water layer. 
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VII. APPENDIX 
A. Fourier transform of Bessel functions 



If 1 ■ 1 /■+°° q\\ 

d^W „2 ^ „2 e"'"" = ^ / ^g|| ^2^„2 -^o(g|ni) 



1 



2^^o(g.|n|) (47) 
1 



(2^) J f„2 , „2A 2^ Jo f„2 , „2 



1 1 ifi(g,||r||) 



(48) 



B. Parameters for static correlation functions 

1. Single bilayer 

For a single supported bilayer and g2|i are the roots of the quadratic equation: 

X X H = 0. 

Solving this equation numerically has to be done with care because of numerical instabilities and we follow the method 
of [48,] and set: 



The roots of the equation are then given by: 



in = — ' 

KlS 



and the coefficients \i , rji are given by: 

2 

Ai = — 3 and A2 = — Ai (50) 

'^111 "'^211) 

771 = 2^2 '72 = -"Hi- (51) 

- 92|| 

2. Double bilayer 

For a double supported bilayer we need to solve a quartic equation leading to q,n (i=1...4). Again, we use the 
method described in [i^. The static coefficients Vi and 77^ ^ are given by: 



(^L-.f,i)(^L-.||,)(^L-.|i,)(^i,2-.|,i) 



(52) 



2 _„2 
2 



"n^d - TT / 2 _ „2 ^ ' (53) 



where Si^i and ^1,2 (resp. ^2,1 and ^2,2) are the solutions of 



K2 A1K2 



and 



The expression for and Ai,j are more complex and will not be given here. 

C. Diagonalization of the free energy for a double supported bilayer 

We have: 

i(ai(g||) + B) -B/2 \fui{q\\ 



The diagonalization of the matrix gives the eigenvalues: 

^+(911 ) = J(ai(9||) + 52(g||) + 2B+^A(9||)) 

= \{^d(l\\) + a2{q\\)+2B- ^ 



and the eigenvectors: 



with: 



In the new basis, 



with: 



TT - B I - 02(911) + A 

" ^ I, ^ ' V 

jj _ B ( ai(g||) -a2(g||) - VA(g||) \ 
A(g||) = (ai(g||)-a2(g||))' + 4B2 

= i^4B2 + - a2(g||) - yA(«^)'. 

^911 = K{i\\)9^{i\\) + ^-{q\\)h'^{i\\), 



^(^11) = 2$(^^ 

- (ai (ay) - 02(911) - yA(^) Ml (gy) +2Bu2(g||) 

'^'"^ = m^) 



We can now apply the theorem of equipartition of energy, 



(iM9ii)n=2A 



2A+ 



21 



ui{q\\) and 'S2('Z||) can now be calculated by inverting Eq. ((58l 



uiiq\\) = ksT- 



U2iq\\) = ksT 



- a2(g||) + v/A(g||)) - {ai{q\\) - ~a2{q\\) - ^ ^{q\\)) ^(g||)g(g||) 



(60) 



D. Coefficients of tiie tliermal correlation functions 



1. Single bilayer 

In the single bilayer case, coefficients of the thermal correlation functions are simply given by: 

1 1 



ai 



9211 



9i|| 



, ^2 



9211 9i|| 



(61) 



2. Double bilayer 

In the double bilayer case, coefficients of the thermal correlation functions are given by: 



a2j 



where /3| ■ are the roots of equations: 



2,2 



]\\> 



(/^M-'?lll)(/^?. 



2-9, 



(62) 

(63) 
(64) 



2 7^ , ^»+^ „ 

X X H — 0. 



(65) 



E. Integration of Eq. (|46|) for off-specular scattering 

When numerically integrating Eq. (|46p cos((j'a;o2;) oscillates very fast whereas f {x,qz) decays logarithmically. Setting 
u — qxgX and next u = v + 2k-K with v S [0; 7r/2] we obtain: 



/o 4^AMr2|rf |i-|'(e,„.e,,)' 1 



Joff-s 



N 

E 



k=0 



Tr/2 



COS (v) e ^'"=0 / 



V + 2fc7r 

92;o 



+ dv COS I w + — I e 

^37r/2 Ai;^''(f + 2fc,r + 7r)'i 

COS {v + tt) e ^'="0 



^ , t; + 2fc7r + 7r/2 
/ 

9^0 



ti + 2fc7r + TT 
/ ,9^ 

9x0 



cos V + 



37r/2 



iV = 100 ensures fast convergence. 



Ag^-^(i. + 2fcir + 37r/2)'' 



, , v + 2kTr + 37r/2 
/ ,9z 

9X0 



(66) 
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F. Structural parameters 





bare 


single 


double 


double 


OTS mixed 




substrate 


bilaycr (20° C) bilayor (49° C) bilayer (62° C) bilayor (42.9° C) 




0.692 


0.692 


0.692 


0.692 


0.692 


crsi[A] 


2 ±0.2 


1.9 ±0.2 


4.6 ±0.2 


3.9 ±0.2 


4 ±0.2 




0.52 ±0.05 


0.52 ±0.05 


0.65 ±0.05 


0.73 ± 0.05 


0.66 ±0.02 




3.8 ±0.5 


3.8 ±0.5 


9.3 ±0.5 


5.1 ±0.5 


6.5 ±0.5 


CSi02 [A] 


2 ±0.2 


1.9 ±0.2 


3.6 ±0.2 


4.8 ±0.2 


4 ±0.2 




1±0.5 


1±0.5 


6.8 ±0.5 


2.8 ± .5 


1±0.1 


Hs 


0.5 ±0.05 


0.5 ±0.05 


0.43 ± 0.05 


0.45 ± 0.05 


0.38 ± 0.05 


Dl,H20[^\ 




7±1 


5.5 ±1 


1.5 ±1 


0.5 ±1 


Pl,head[e~ /A'^] 




0.45 ± 0.02 


0.42 ±0.02 


0.40 ± 0.02 


0.5 ± 0.02 


pi.taule" /A'^] 




0.33 ± 0.02 


0.30 ± 0.02 


0.32 ±0.02 


0.31 ±0.02 






0.26 ± 0.02 


0.28 ±0.02 


0.25 ± 0.02 


0.26 ± 0.02 


P2,tail[e'/A^] 




0.33 ±0.02 


0.29 ±0.02 


0.30 ± 0.02 


0.29 ± 0.02 


P2,head[e^ / A^] 




0.44 ± 0.02 


0.41 ±0.02 


0.46 ±0.02 


0.51 ±0.02 


dl^head[A] 




3.8 ±0.5 


6.1 ±0.5 


5.1 ±0.5 


7 ±0.5 


Dl,tail[A] 




22.5 ± 0.5 


26.0 ± 0.5 


26.9 ± 0.5 


22.7 ±0.5 


dcH3,Bi [A] 




2.0 ±0.5 


1.5 ±0.5 


1.75 ± 0.5 


2.1 ±0.5 


D2,taillA] 




23.0 ± 0.5 


24.5 ± 0.5 


25.7 ±0.5 


21 ± 0.5 


d2,head[A] 




4.3 ±0.5 


6.1 ±0.5 


5.7 ±0.5 


3.5 ±0.5 








28.0 ± 1 


25.6 ± 1 


25.1 ± 1 


P3,hea.d[e~ / A^] 






0.41 ±0.02 


0.40 ±0.02 


0.39 ±0.02 


pZ,tail [e"M^] 






0.29 ± 0.02 


0.25 ± 0.02 


0.28 ±0.02 


PCH3,B2 [e~/A^] 






0.26 ± 0.02 


0.24 ± 0.02 


0.26 ± 0.02 


PA,tail[e-~ 1 






0.31 ± 0.02 


0.33 ± 0.02 


0.32 ± 0.02 


pA,h^ad{e~ l 'A^^\ 






0.42 ±0.02 


0.41 ±0.02 


0.40 ± 0.02 








3.1 ±0.5 


3.5 ±0.5 


2.9 ±0.5 








20.7 ±0.5 


16.3 ±0.5 


21.0 ±0.5 


dcH3,B2 [A] 






2.1 ±0.5 


2.2 ±0.5 


1.8 ±0.5 


D4,tail [A] 






20.6 ± 0.5 


19.6 ±0.5 


21.6 ±0.5 


d4,head[A] 






3.4 ±0.5 


3.9 ±0.5 


3.9 ±0.5 


71 [mN/m] 




70 ±10 


15.0 ± 1 


5.0 ±0.2 


15.8 ±2 






350 ± 100 


370 ± 50 


110 ±25 


350 ± 100 


Logim%[J.m-^]) 




11.2 ± 1 


10.3 ± 1 


10.9 ± 1 


10.3 ± 1 


72 [mN/m] 






1.0 ±1 


0.09 ± 0.1 


0.1 ±0.05 


K2[kBT] 






150 ± 50 


45 ± 25 


200 ± 50 


Log{m%[J.m-^]) 






8.9 ± 1 


8.3 ±1 


10.0 ± 1 


Log{m',2[J-m-^]) 






12.2 ±0.3 


13.4 ± 1 


12.2 ± 1 



TABLE I: p: electron density; D: box thickness; d: Gaussian width; a: roughness; ^sicutoff, Hg: roughness exponent, 7: 
tension, k: bending rigidity, U": second derivative of the interaction potential. Subscript s represents substrate, 1 membrane 
1 (next to the substrate), and 2 membrane 2. 
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